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ABSTRACT 

To understand the physics of solar flares, including the local reorganisation 
of the magnetic fleld and the acceleration of energetic particles, we have flrst to 
estimate the free magnetic energy available for such phenomena, which can be 
converted into kinetic and thermal energy. The free magnetic energy is the excess 
energy of a magnetic conflguration compared to the minimum-energy state, which 
is a linear force- free fleld if the magnetic helicity of the conflguration is conserved. 
We investigate the values of the free magnetic energy estimated from either the 
excess energy in extrapolated flelds or the magnetic virial theorem. For four 
different active regions, we have reconstructed the nonlinear force-free field and 
the linear force-free fleld corresponding to the minimum-energy state. The free 
magnetic energies are then computed. From the energy budget and the observed 
magnetic activity in the active region, we conclude that the free energy above the 
minimum-energy state gives a better estimate and more insights into the flare 
process than the free energy above the potential fleld state. 

Subject headings: Sun: magnetic flelds — Sun: flares — Sun: magnetic energy 
— Sun: magnetic helicity 



1. Introduction 

Due to the low value of the plasma (3 (the ratio of gas pressure to magnetic pressure), the 
solar corona is magnetically dominated. To describe the equilibrium structure of the coronal 
magnetic fleld when gravity is negligible, the force-free assumption is then appropriate: 

VAB = aB, (1) 

where a = gives the potential (or current-free) fleld, a = est gives the linear force-free fleld 
[Iff), and a being a function of space gives the noi ilinear force- free fleld (nlfD . The properties 



of force-free flelds have been well described (e.g., IWoltjerlll958l : lMolodenskiilll969 



Alvl 11984 
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Bergerlll985l ). IWoltjerl (119581 ) with general astrophysical configurations in mind derived two 
important theorems: (i) in the ideal MHD limit the magnetic helicity is invariant during the 
evolution of any closed flux systems, (ii) the minimuna energy state i s the li r iear f orce-free 
field conserving the magnetic helicity (see also lAlylll984j : lBergerlll985l ) . iTaylorl (119861 ) applied 
this to laboratory experiments and hypothesized that in a weak but finite resistive regime 
the total magnetic helicity of the flux s ystem is i nvaria nt during the relaxation process to 
a minimum energy state. According to IWoltjen (119581 ). the relaxed state is then a linear 
force-free field. Therefore the free magnetic energy that can be released during a relaxation 
process is the excess energy of the magnetic configuration above the linear force-free field 
with the same magnetic helicity. 



Heyvaerts &: PriestI (119841 ) were the first to suggest the importance of magnetic helicity 
and Taylor relaxation in the solar corona. They extended the Woltjer- Taylor theory for an 
isolated structure bounded by magnetic surfaces to one that of a coronal field in which the 
field lines enter or leave the volume (through the photosphere): thus the magnetic helic- 
ity is allowed to enter or leave the corona as the photospheric field changes in time. They 
also suggested that the coronal field evolves locally through a set of linear force-free fields 
with the field continually relaxing and the footpoint connections continually changing by 
small-scale turbulent reconnections, which heat the corona. Moreover they suggested that, 
if the magnetic helicity becomes too large, an eruption takes place in order to expel the 
excess magnetic helicity. The coronal heating mec hanism by magnetic turbu lent relaxation 
was later developed into a self-consistent t heorv (jHevvaerts fc PriestI Il992l ). Based on a 
statistical analysis of vector magnetograms, iNandy et al.l (120031 ) have shown that the relax- 
ation process of fia r e-prod uctive active regions is similar to Taylor's theory. Nevertheless, 
Regnier fc Canfieldl (120061 ) have shown that the magnetic helicity can evolve significantly on 
a short time scale (about 15 min) and that the evolution of the coronal magnetic field is 
often well described by a series of nonlinear force-free equilibria. The modelled evolution 
of global coronal fields by s u ccessiv e nonlinear force-free equilibria was also investigated by 
Mackay fc van Ballegooijen ( 2007al lbl). 



To better understand the physics of flares, we need to estimate the amount of mag- 
netic energy available in a magnetic configuration for conversion into kinetic energy and/or 
thermal energy in a solar flare. There is no free magnetic energy in a potential field configu- 
ration: this is a minimum-energy state for a given normal magnetic field at the photosphere, 
and the magnetic energy depends only on the distribution and amount of flux through the 
photosphere. The linear and nonlinea r configurations, howeve r, do have free energy due to 
the presence of currents. As shown in lRegnier fc PriestI (120071 ). the energy storage in active 
regions can be (i) in the corona due to the existence of large-scale twisted flux bundles, or 
(ii) near the base of the corona associated with the existence of a complex topology. The 
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free energy can be estimated from photospheric orchromospheric magnetic fields based on 
the magnetic virial theorem (IMolodenskiil Il969l : lAlylll984j ). or from reconstructed 3D coro- 
nal fields (often assuming a force-free equilibrium). Using nonlinear force-free fields, the 



magnetic energy budget has been estimated before and after a flare (iBleybel et al.l 12002 



Regnier fc Canfieldl 120061 ): as expected the authors found that the magnetic energy usually 
decreases during the flare. Nevertheless, it strongly depends on the strength of the flare, on 
the processes of energy injection (e.g., flux emergence, flux cancellat i on, su nspot rotation) 
and on the time span between the reconstructed fields. iBleybel et al.l (120021 ) have suggested 
that Taylor's theory does not apply to flares and C MEs. The same conclu sion has been 
reached previously by numerical simulations (see e.g.. lAmari fc Lucianil |2000| ) . This can be 
understood if (i) the helicity is not conserved during a flare or a CME in the finite domain 
of computation due to the injection of helicity through the photosphere or into the CMEs, 
(ii) the eruption phenomenon is often localized in the active region and so does not affect or 
modify strongly the nonpotentiality of the field outside the flare surroundings. Note that the 
energy flux (or Poynting flux) can be deri ved from successive magnetic field measurements 
when the plasma flows are known (see e.g.. lKusano et al.ll2002l ). The Poynting flux gives an 
estimate of the injected energy through the photospheric surface due to transverse motions 
and/or flux emergence. 

In this letter, we compute the free magnetic energy for different active regions assuming 
a nonlinear force-free equilibrium with a reference field being either the potential field or the 
linear force-free field, and from the magnetic virial theorem. We are assuming that at a give n 
time and with the same boundary conditions the minimum-energy state is given (jWoltjerl 's 
theorem) by the linear force-free field with the same magnetic helicity as the nonlinear force- 
free field. We are not here investigating the validity of the Taylor- Hey vaerts theory in solar 
active regions. 



2. Selected Active Regions 

In order to compare the different measurements of free magnetic energy, we have selected 
four different active regions with different types of activity (confined flares, flares associated 
with a CME or filament eruptions) and at a different stage of their evolution (before or after 
a flare): 



AR8151: observed on February 11, 1998 at 17:36 UT, this is an old decaying active region 
(decreasing magnetic flux and magnetic polarities diffusing away). A filament eruption 
associated with an aborted CME was reported on Feb. 12, but no flare was observed. The 
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vecto r magnetic field was recorded by the MEES/IVM (IMickey et al.lll996l : iLaBonte et al. 



19991) ■ The high values of the current density imp ly strongly sheared and twisted flux bundles 



see 



Regnier et al.ll2002l : iRegnier fc Amarill2004l ). Due to the existence of highly twisted flux 



tubes (with more than 1 turn) and the stability of the reconstructed filament and sigmoid 
(with less than 1 turn), the authors concluded that the eruptive phenomena was most likely 
to be due to the development of a kink instability in the highly twisted flux bundles; 



AR8210: observed on May 1, 1998 from 17:0 to 21:30 UT, this is a ne wly emerged active 
region with a complex topology as described in [Regnier fc Canfieldl fl200fih . A Ml. 2 flare was 
recorded on May 1, 1998 at 22:30 UT. The selected vector magnetogram ( MEES/IVM) at 
19:40 UT was observed during a "quiet" period between two C-class flares. In lRegnier fc Canfield 
( 120061 ). the authors described the magnetic reconnection processes occurring during this time 
period and leading to a local reorganisation of the magnetic field. The reconnection pro- 
cesses are related to the slow clockwise rotation of the main sunspot or a fast moving, newly 
emerged polarity. Following the time evolution during 4 hours, the authors showed that the 
free magnetic energy decreases during the flare over a period of about 15 min, and the total 
magnetic energy is slightly increased during this time period; 



AR9077: this corresponds to the famous Bastille day flare in 2000 ( 


Liu & Zhane 


2001 


Yan et al. 


2001; 


Fletcher & Hudson 


2001; 


Kosovichev & Zharkova 


2001; 


Wang et al. 


2005) 



The vector magnetogram was recorded at 16:33 UT after the X5.7 flare which occurred at 
10:30 UT. The active region was still in the magnetic reorganisation phase after the flare and 
"post" -flare loops were observed in 195A TRACE EUV images. The flare was also associated 
with a CME; 



AR10486: this active region is responsible for the main eruptions observed during the 
Halloween events (26 Oct. to 4 Nov. 2003). The MEES/IVM vector magnetogram was 
recorded on October 27, 2003 at 18:36 UT before the X17.2 flare which occurred at 11:10 
UT on October 28. The flaring activity o f this active region and the associated CMEs 
have been extensively studied. For instance, iMetcalf et al.l (120051 ) have shown that the large 
magnetic energy budget (~ 3 10'^^ erg) on Oct. 29 is enough to power the extreme activity 
of this active region. 

For these particular active regions, the reduction of the full Sto kes vector to derive 
the magnetic field has a lready been detailed in several articles (e.g., iRegnier et al.l 12002 



Regnier fc Canfieldll2006l ) - the 180- degree ambiguit y in th e transver se component was so lved 
by using the algorithm developed in lCanfield et al.l (119931 ) (see also lMetcalf et al.ll2006l ). 
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3. Magnetic Fields 

From the observed vector magnetic field as described in Section [21 we extrapolate to 
obtain three types of coronal magnetic field, each of which has the vertical component of the 
magnetic field imposed at the photopshere: 

- potential field: there is no current fiowing in the magnetic configuration; this is the 
minimum-energy state that the magnetic field can reach when the magnetic helicity is 
not conserved; 

- linear force-free field: we compute the linear force-free whose a parameter is chosen so 
that the total magnetic helicity is the same as the nonlinear force-free; in other words, 
this gives the minimum-energy state that conserves the magnetic helicity; 



nonlinear force- free field : we use the vector potential Grad- Rubin-like method (iGrad fc Rubin 



19581 : lAmari et al.lll999l ). The bottom boundary conditions also require the knowledge 



of a in one polarity derived from the transverse field components: a = ^ ~ ■ 



In order to have energy values which can be compared, we have imposed the same 
closed conditions on the side and top boundaries for each model. To satisfy these conditions, 
we surround the vector magnetic field observed by MEES /IVM by weak field measurements 
provided by SOHO/MDI line-of-sight observations. The active region fields are then confined 
by a surrounding potential field and the magnetic field decreases from the center of the active 
region (compatible with the field vanishing at infinity). The magnetic fiux is balanced. 



4. Free Magnetic Energy 

From the 3D coronal magnetic configurations, we can derive the magnetic energy for 
the different active regions and different models: 

Em= I (2) 

in a volume Q. The free magnetic energy is derived from the nonlinear [nlff) force-free field 
using either the potential (pot) or linear force-free (Iff) field as reference field: 

A</^ = - (3) 
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and 



(4) 



The Iff field used here has the same relative magnetic helicity as the nlff field satisfying 
Woltjer's theorem. That implies that the nljf field has to be computed first, and then the 
Ijf field is determined by an iterative scheme to find the a value matching the hel i city o f 
the nlff field. Th e relative magnetic helicity is computed from the iBerger &: Field! (119841 ) 
equation (see e.g. 



Regnier et all 120051 ): 



/ (A-Apot) 
Jn 



[B + Bpot) dn 



(5) 



where B and A (resp. Bpot and Apot) are the nlff (resp. potential) magnetic field and its 
associated vector potential computed in the volume Q. The relative magnetic helicity given 
by Eqn. ([5]) satisfies the closed boundary conditions used by the Grad-Rubin reconstruction 
method. 



Table 1: Free magnetic energy and relative magnetic helicity for the different active regions 





m 

(1032 erg) 


(10^2 erg) 


a 

(Mm-i) 


A /^"'^'-^ 
Iff 

(10-32 erg) 


AEZ 1 
observed 


;i032 erg) 
computed 


AHm 
(10^2 Mx2) 


Comments 




AR8151 


0.64 


0.26 


0.067 


0.05 


1.2 


0.79 


0.47 


twisted bundles 


-a 


AR8210 


10.6 


0.24 


-0.056 


0.14 


1.63 


0.79 


-4.2 


before C flare 


1 


AR9077 


14.2 


2.21 


-0.015 


1.62 


0.48 


1.25 


-14.6 


post-flare loops 




AR10486 


70.5 


18.05 


0.021 


7.23 


41.7 


2.62 


35.1 


before X17 flare 
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For the sake of comparison, we also compute the free niagnetic energy derived from the 



magnetic virial theorem assuming a force-free field ( e. g. lAlyl Il989l : iKlimchuk et al.l Il992 



Metcalf et al.lll995l . l2005l : IWheatland fc Metcalj 120061 ). Considering that the magnetic field 

Spot + b, then following 



can be dec omposed into a potential part and a nonpotential one, B 
Alyl (119891 ) the free magnetic energy (above potential) is: 

1 

' 47r 



AEi: 



{xhx + yhy)B^ dxdy 



(6) 



in the half- space above the surface S. The free magnetic energy from the virial theorem only 
requires the magnetic field distribution on the bottom boundary. We compute Eqn. from 
either the observed vector magnetic field (not necessarily force-free) or the reconstructed 
nlff field on the photosphere. It is important to note that the energy values derived from 
th e magnetic viria l theor em are strongly influenced by the spatial resolution as mentioned 



m 



Klimchuk et al.l fll992h . 



Considering that the Iff field is the minimum energy state of the n{ 
values can be sorted as follows: 



'field, the energy 



(7) 



In Fig. [H we plot the free energy values in the reconstructed magnetic configurations 
using the potential field as reference field for the nlff fields (triangles) and Iff fields (crosses). 
The difference between the two values is the minimum free energy AE^jj-^ according to 
Woltjer's theorem. Figure [1] clearly shows that the free magnetic energy can vary by at 
least 2 orders of magnitude: the energy is strongly influenced by the total magnetic flux and 
the distribution of the polarities. By comparing the amount of free energy A£;;^f^ and the 
observed eruptive phenomena, we can conclude that AE^^j-^ gives a better estimate of the 
free energy. For instance, AEpg{^ is similar for AR8151 and AR8210 but AE^jj-^ is nearly 
three times larger for AR8210. And the related eruptive phenomena are very different: a slow 
fllament eruption without a flare for AR8151 and a C-class flare for AR8210. For AR9077, 
AE"jj-^ is still enough to trigger an X-class flare but certainly not the X5.7 flare observed 
prior to the time considered here. For AR10486, AE^^j-^ is signiflcantly reduced compared 
to AEpjj{-^ but still e nough to trigger pow erful flares which explains the high level of activity 



in this active region (IMetcalf et al.ll2005l ). 



In Table [T], we summarize the different values of free magnetic energy, the magnetic en- 
ergy of the n/j^ magnetic conflgurations (i?^'-^-^) and the relative magnetic helicity. We also 
mention the a values used to compute the Iff fields satisfying Woltjer's theorem. We notice 
that the different values of free energy are consistent and increase when the eruption phenom- 
ena increase in strength with the exception of AE^^ from the observed magnetograms. The 
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Fig. 1. — Magnetic energy above potential for both the nlff Held (triangles) and the Iff Held 
(crosses) of the four selected active regions (units of 10^^ erg). The free magnetic energies 
AE^^j^ are given by the differences between the triangles and crosses. 
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latter is related to the applicability of the vi rial theorem because the observed magnetograms 
are not force- free at the photospheric level (iMetcalf et al.lll995l ). 



Conclusions 



We have computed the free magnetic energy from several formulae in various active 
regions at different stages of their evolution: from the difference between the nljf field and 
either the potential field {AE^jjf-^) or the Iff field {AE^jj-^) having the same magnetic helicity, 
and from the magnetic virial theorem {AE^"^) using either the observed field or the n/jf field. 

The free magnetic energy AE^jj^ is a better estimate (than AE^H^) of the energy budget 
of an active region available for flaring assuming that the magnetic helicity is conserved 
and gives more insights into the possible eruption mechanisms in the active region. For 
AR8151, it is clear that there is not enough energy to trigger a flare c apable of a large 



scale reorganisation of the field (~ 5 lO^'' erg). Therefore as stated in iRegnier fc Amari 



(120041 ) the kink instability of the highly twisted flux tube is most likely to be responsible 
for the observed eruptive phenomenon. Despite a magnetic energy of about 10^^ erg, the 
free magnetic energy in AR8210 is only 1% of the total energy but is enough to trigger 
small confined flare s . Thi s is consistent with the observations and modelling described in 



Regnier fc Canfie. 



(iRegnier fc Priest 



dl (120061) . We note that for the two possible mechanisms to store energy 



20071 ). the presence of large twisted flux bundles is more efficient than the 



highly complex topology: 10% of free energy in AR8151 compared to 1% for AR8210. The 
magnetic energy budget of AR9077 is still important even if the observed field is after a X5.7 
flare. Therefore even after a strong flare with post-flare loops resembling potential field lines, 
the magnetic configuration is far from potential and the energy budget is still sufficient to 
trigger further powerful flares. For AR10486, AE^^j^ is certainly not sufficient to trigger the 
observed XI 7. 2 flare, but the AEpjj{^ seems to be more consistent with the recorded flaring 
activity. This can be explained by the fact that the main hypothesis of Woltjer's theorem 
is not satisfied: the X-class flare is associated with a CME expelling a magnetic cloud (and 
therefore magnetic helicity) into the interplanetary medium. 

The free magnetic energy AE^^^ gives consistent values when computed from the n/jf ex- 
trapolated fields. For most photospheric magnetograms, the force-free assu mption is not 
well satisfied and so leads to inaccurate values of AE^!^"^ 



from observations. In Metcalf et al. 



(120051 ) ■ the computation of the free energy from the virial theorem was performed using 
chromospheric magnetic field measurements which a re more force- free than photospheric 
magnetograms (IMetcalf et al.lll995l : iMoon et al.ll2002l ). 
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To have a better understanding of flaring activity, our main conclusion is that it is useful 
to compute both AE'^j^f-^ and AE^jj-^: the first giving an upper limit on the magnetic energy 
that can be released during a large flare, especially when associated with a CME, the second 
being a good estimate of the energy budget for small flares and allowing us to distinguish 
between different flare scenarios. 

We thank the UK STFC for financial support (STFC RG). The computations were 
done with XTRAPOL code developed by T. Amari (supported by the Ecole Polytechnique, 
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